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ABSTRACT 

Recent studies of eukaryotic DNA replication timing 
profiles suggest that the time-dependent rate of 
origin firing, /(f), has a universal shape, which 
ensures a reproducible replication completion 
time. However, measurements of /(f) are based on 
population averages, which may bias the shape of 
the /(f) because of imperfect cell synchrony and 
cell-to-cell variability. Here, we measure the 
population-averaged /(f) profile from synchronized 
Saccharomyces cerevisiae cells using DNA 
combing and we extract the single-cell /(f) profile 
using numerical deconvolution. The single cell /(f) 
and the population-averaged /(f) extracted from 
DNA combing and replication timing profiles are 
similar, indicating a genome scale invariance of 
the replication process, and excluding cell-to-cell 
variability in replication time as an explanation for 
the shape of /(f). The single cell /(f) correlates with 
fork density in wild-type cells, which is specifically 
loosened in late S phase in the clb5A mutant. A pre- 
viously proposed numerical model that reproduces 
the wild-type /(f) profile, could also describe the 
clb5A mutant /(f) once modified to incorporate the 
decline in CDK activity and the looser dependency 
of initiation on fork density in the absence of Clb5p. 
Overall, these results suggest that the replication 
forks emanating from early fired origins facilitate 
origin firing in later-replicating regions. 

INTRODUCTION 

Eukaryotic DNA replication starts at multiple origins that 
are fired at different times during S phase (1). The mech- 
anisms that control the temporal order of origin firing 
through S phase are not clearly understood. A topic of 
recent debate is whether origins fire independently of 



each other with variable efficiencies that define their 
mean activation time (2,3), or are coordinated by some 
mechanism that subordinates the firing of late origins to 
the progression of DNA synthesis in earlier replicons 
(4,5). Czajkowsky et al. (6) visualized by DNA combing 
the replication patterns of 5. cerevisae chromosome VI. 
Contrary to what would be expected for a strictly 
regulated temporal program, no two molecules exhibited 
the same pattern and replication of different regions of the 
same chromosome occurred independently of each other. 
Nevertheless, averaging the patterns of all examined fibres 
recapitulated the population-averaged replication timing 
profiles obtained previously (7-9). These results favour a 
probabilistic model for replication timing control. On the 
other hand, McCune et al. (10) showed that in the clb5A 
mutant, which affects the origin firing activity of 
Cyclin-Dependent Kinase (CDK) late in S phase, the rep- 
lication timing profile is only affected in the regions that 
replicate late in the wild-type, which replicate even later in 
the mutant. Because the efficiency of early origins was not 
affected, McCune et al. disfavoured a purely stochastic 
model of origin firing. 

An essential kinetic parameter of S phase is the time- 
dependent rate of firing, I(t). Single-molecule analysis 
of replication intermediates in Xenopus egg extracts 
have shown that /(/) increases during S phase to allow 
the timely completion of DNA replication (11-15). 
Subsequent work using both DNA combing in Xenopus 
and analysis of whole genome replication timing profiles 
in a wide range of other eukaryotes confirmed that /(/) 
increases during S phase but also suggested that I(t) de- 
creases before the end of S phase in a similar manner in all 
these organisms (16). However, was extracted from 
replication timing profiles by counting peaks whose 
height was read as the origins firing time. As these 
profiles are averaged over large cell populations, which 
masks cell-to-cell variability, the height of peaks in fact 
depends on both origin firing time and origin firing effi- 
ciency (17,18). Thus, counting visible peaks as origins 
leads to overestimate both their firing rate and firing 
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time. On the other hand, inefficient origins fail to trace 
detectable peaks. If later origins are less efficient, this may 
lead to underestimate the rate of origin firing late in S 
phase. A more sophisticated mathematical analysis of a 
Saccharomyces cerevisiae genome-wide replication time 
course profile indeed suggested that the decrease in /(/) 
was less prominent than initially measured (3). Finally, 
measurements based on population averages may particu- 
larly affect the aspect of the I(t) profile at the start and the 
end of S phase because of imperfect cell synchrony. 

Assuming that the population-averaged /(/) is a reason- 
able approximation of the single cell I(t), two models have 
been proposed to quantitatively account for the reported 
profiles (19,20) (reviewed in 2 and 21). Both the models 
require that the amount of a replication factor increases 
during S phase to fully explain the increasing part of the 
profile (a constant amount of the replication factor 
would also result in an increasing I(t), but the experimen- 
tal profile increases even more strongly, which the authors 
attribute to an increase in the replication factor). They 
differ in their explanation of the decreasing part. Goldar 
et al. (20) suggest that it results from a coupling between 
the density of replication forks and the probability of 
origin firing at unfired origins. Gauthier and Bechhoefer 
(19) suggest that the decreasing part of I(t) can be ex- 
plained if the search process between DNA and some ini- 
tiation factor is subdiffusive, without explicitly tying firing 
rate to fork density. 

Transacting factors that control replication initiation 
include CDK and Dbf4-dependent kinases (22). The 
effect of altering CDK activity during S phase in S. 
cerevisiae has been analysed in detail (23-26). During a 
normal S phase, CDK is activated by Clb5p or a 
second, short-lived, cyclin Clb6p, to stimulate origin 
firing. The absence of Clb5p delays the replication of 
multiple zones that together compose over half of the 
genome, due to decreased origin firing and elicits activa- 
tion of a Rad53-dependent checkpoint response that sta- 
bilizes replication forks (10,27). Strikingly, the delayed 
replication zones in the clb5A mutant coincide with 
delayed replication zones seen during replication stress 
in a mutant lacking the Isw2 and Ino80 chromatin 
remodelling complexes, which are normally enriched at 
replication sites (28). These results suggest that chromatin 
remodelling activities located at replication forks specific- 
ally promote fork progression and/or origin firing in those 
late-replicating regions. 

In this work, we deliberately compare in a wild-type 
S. cerevisiae strain the I(t) obtained from combing data, 
where errors due to variable origin efficiency do not 
exist, with an /(/) derived from replication timing 
profiles as in previous works. We then extract from the 
combing data the average single-cell /(/) profile. We show 
that the three profiles are very similar and that the firing 
rate correlates with the fork density in wild-type cells. This 
dependency is specifically loosened in late S phase in the 
clb5A mutant. The model previously proposed by Goldar 
et al. (20), which accounts for the wild-type /(/) profile, 
could also describe the clb5A mutant I(t) profile once the 
decline in CDK activity and the looser dependency of 
origin firing on fork density in the absence of Clb5p 



were incorporated into the model. Overall, these results 
suggest that replication forks emanating from early- 
activated origins facilitate origin firing in later-replicating 
regions. 

MATERIALS AND METHODS 

The strain 

The strain used in this study, MCM869 (MATa ade2-l 
trpl-1 canl-100 leu2-3, his3-ll,15 URA3::GPD- 
TK7x AuRlc::ADH— hENTl barlA::LEU2 cdc21A:: 
kanMX), is a derivative of the E1000 strain (29) in 
which the CDC21 gene (encoding thymidylate synthase) 
was disrupted to increase BrdU incorporation (30). 
During growth, MCM869 is maintained in medium 
containing 100 uM thymidine. Unless indicated other- 
wise, cells were grown at 30°C in synthetic minimal 
media (31) supplemented with adenine and casamino 
acids (2%). 

Cell cycle analysis 

Approximately, 2 x 10 9 cells were arrested in Gl phase by 
incubation for 4 h in the presence of 1 |iM oc-factor and 
100 uM thymidine (Sigma), centrifuged and washed twice 
with fresh medium, transferred to fresh medium contain- 
ing 1 uM a-factor and lOOuMBrdU for 30min, released 
into S phase by washing out the a-factor by centrifugation 
and transferring in fresh medium containing 100 uM Brdu 
and 30ug/ml nocodazole (Sigma). Samples (2 x 10 8 cells) 
were collected every 15min over a 150-min period, at 
which time the cells reached the G2/M phase. An 
aliquot (2 x 10 4 cells) of each timed sample was removed 
to monitor S phase progression by flow cytometry accord- 
ing to Epstein and Cross (23) and the rest was centrifuged 
and washed twice and resuspended in fresh medium con- 
taining 30ug/ml nocodazole and 100 uM thymidine at 
30°C until the G2/M phase. Each sample was divided in 
two parts, one (15 x 10 7 cells) was analysed by DNA 
combing and the other (5 x 10 7 cells) was used to 
quantify BrdU incorporation. 

Quantitation of BrdU incorporation 

Genomic DNA was prepared as described by Sambrook 
and Russell (32). For each time point, equal amounts of 
genomic DNA (100 ng) were denatured in 0.4 M NaOH 
(Sigma-Aldrich), spotted onto a nitrocellulose membrane 
(GE Healthcare) and UV-cross linked. Immunodetection 
was performed using a monoclonal anti-BrdU antibody 
(Becton-Dickinson) and a secondary IgG coupled to 
horseradish peroxidase (Ebioscience). Signals were 
detected by electrochemiluminescence and quantified 
with ImageJ software. Serial dilutions of BrdU-labelled 
DNA were systematically performed to check the linearity 
of the signal. The copy number of DNA y was deduced 
from the fraction of incorporated BrdU, .\- Brdu 
(normalized by the maximal amount of incorporated 
BrdU) as y = 1 +x Brdu . The fraction of replicated DNA 
/ona represents the amount of replicated DNA divided by 
the total amount of DNA present initially. The fraction of 
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incorporated BrdU, XBrdu> corresponds to the amount of 
BrdU incorporated in to the double stranded DNA 
divided by the total amount of DNA. Assuming that the 
length of the DNA is L, then the total amount of DNA is 
equal to L + Lx/ DNA . So x Brdu = 2 x/ DNA /(l +/ DNA ) 
and therefore, / DNA = x BrdU /(2 - x Brdu ). 

DNA combing 

Silanized glass coverslips were prepared as described (33) 
and used to comb BrdU-labelled DNA from the 60, 90 
and 105 min time points as described (34). Total and 
BrdU-labelled DNA was detected with fluorescent 
antibodies as described (34). Total DNA was stained in 
red using an anti-DNA antibody (Abcys) and a secondary 
IgG coupled to Alexafluor 532 (Invitrogen). BrdU tracks 
were stained in green using a monoclonal anti-BrdU 
antibody (Becton-Dickinson) and two secondary IgGs 
coupled with Alexafluor 488 (Invitrogen). Labelled fibres 
were imaged using a Leica epifluorescence microscope 
equipped with a x 100 objective and a charge-coupled 
device camera. The fluorescence signals were analysed 
using a Matlab-(Mathworks) based homemade software. 
Green tracks <6 kb were counted as recently fired origins 
(tracks <0.5kb were excluded to avoid false positives due 
to background fluorescence). To construct the I(t) profile, 
it would seem straightforward to use the BrdU chase time. 
Unfortunately, there is a broad dispersion of S phase 
starting time in the cell population, as shown by the 
flow cytometry data, so it is impossible to assign to each 
fibre the exact time at which the cell it comes from, started 
replicating. On the other hand, the local time at which 
each fibre started replicating can be inferred relatively pre- 
cisely using /dna(0- Therefore, the combed DNA fibres 
were sorted according to replication extent / DNA . For each 
fibre, eyes 0.5-6.0 kb long (detectable initiation events that 
occurred over a 2-min interval, assuming that the speed of 
replication fork is 1.46kb min -1 (35)) were counted and 
the number was divided by 2 min and by the length of 
unreplicated DNA, to derive the frequency of origin 
firing as a function of / DNA . The 7(/dna) data points 
were converted to I(t) using the /dna(0 curve determined 
from kinetic measurements of replicating cells. This 
strategy implicitly assumes that the I(t) is spatially homo- 
geneous and was posteriorly confirmed by the collapse of 
the local and global /(/) profiles. The fork density was 
obtained by counting the number of forks and dividing 
by the length of each combed DNA fibre. The complete 
characterization of all fibres analysed in this study is 
described in Supplementary Section 5. 

Quantitative flow cytometry analysis 

To extract the rate of cell DNA synthesis from the distri- 
bution of DNA content over the population we used the 
mathematical framework developed by Bertuzzi et al. (36). 
This analysis is based on two hypotheses: the uniform law 
of cell-DNA synthesis over the population and the 
absence of cell losses. The validity of the first assumption 
is further discussed in Supplementary Section 5. Adapting 
equation (14) of Bertuzzi et al. (36) to our study, we obtain 



for the fluorescence density normalized to unit area, v (t;,f) 
at a time t: 

v&t) = v(m+m, i)[2(i - #0 - h(t)]+h(t) 

J °, (1) 

h(t)= f /(r)^,l+/DNA(r))dr 

Jt-T, 

h{t) = / /(r)dr 
Jo 

where £ is the position of the detection channel 
normalized by the position of the Gl peak, 0 X is the 
fraction of cells in Gl and the function _AT(£,7+/" dna (t)) 
is the kernel function (37) representing the dispersion 
of propidium iodide (PI) fluorescence over the detection 
channels. We then extracted the fraction of cells in 
S phase, n(t), as follows. We first determined the 
population-averaged kinetics of DNA replication, 
/dna(0 an d the total rate of DNA synthesis, 
d/DNA(0/d?- We then extracted from the flow cytometry 
profiles the S phase influx function, fir), fir) is continuous 
and differentiate at each of its points and therefore can 
be considered as a sum of continuous and differentiable 
functions, fit) = 1/6 J^=\ A i /V2n(r i e- ( '-'' )2/2c '" . Therefore, 
fit) can be extracted by fitting all six flow cytometry 
profiles simultaneously using equation (1) and a simplex 
algorithm. The time-dependent distribution of cellular 
DNA content, n(x,t), was then calculated as 
/(r)/(fi/DNAW/dr) (Figure 2C). Finally, the probability 
for a cell to be in S phase is calculated as 
n(t) = /j n(x,t)dx/Ntoub where iV tota i represents the total 
number of cells. 

Analysis of published microarray data 

The time-dependent rate of origin firing I(t) in the budding 
yeast S. cerevisiae was constructed using the genome-wide 
replication profiles determined by Raghuraman et al. (7), 
Yabuki et al. (8) and McCune et al. (10), where the mean 
replication time of each studied locus is plotted against 
chromosomal position so that efficient origins appear as 
peaks. The replication profile of each chromosome was cut 
into 1-min slices of time. I(t) was calculated for each time 
point as the ratio between the number of fired origins and 
the amount of unreplicated DNA. The number of active 
replication forks at time t was detected as the intercept 
between the replication profile curve and a y = t horizon- 
tal line. The replication fork density was calculated as the 
ratio between this number and the genome length. The 
length of S phase of our combing experiments (synchron- 
ization using a-factor, S phase at 30°C; T s = 42 ± 15 min) 
was comparable with that measured by Raghuraman et al. 
(43 min) and McCune et al. (45 min; in both the studies 
double synchronization by a-factor and cdc7 mutant was 
used, S phase at 30°C) but longer than that measured by 
Yabuki et al. (8) (synchronization using a-factor, S phase 
at 26°C; 30 min). When the firing time for known ARSs 
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Figure 1. Experimental workflow and flow cytometry analysis of DNA 
content following a-factor release. (A) Cells were synchronized in Gl 
with a-factor and synchronously released into S phase by washing out 
the a-factor and transferring in fresh medium containing BrdU and 
nocodazole. Samples were collected every 15 min over a 150-min 
period, at which time the cells reached the G2/M phase. An aliquot 
of each timed sample was used to monitor S phase progression by flow 
cytometry fluorescence-activated cell sorter. The remaining part of the 
timed sample was transferred to fresh medium containing nocodazole 
and thymidine at 30°C up to 150 min. Each sample was divided in two 
parts, one was analysed by DNA combing and the other part was used 
to quantify BrdU incorporation. (B) Experimental data (open symbols) 



(from OriDB database (38)) measured by Raghuraman 
et al. was plotted against the firing time measured 
by Yabuki et al., a linear relation was observed 
[? ya buki= 1 3.26 (± 0.68) + 0.43(± 0.03) / Raghuraman , P< 
10 4 ]. Therefore, we used the S phase duration of 
Raghuraman et al. as our reference point and corrected 
the timing measured by Yabuki et al. by the obtained 
linear relation. 



RESULTS 

To study the dynamics of origin firing, Gl -arrested cells 
were synchronously released into S phase in the presence 
of BrdU at / = 0. Every 15 min over a 150min period, 
three samples were harvested from the solution. The first 
sample was analysed by flow cytometry, the second sample 
was used to measure the amount BrdU incorporated into 
the genomic DNA and the third sample was analysed by 
DNA combing (Figure 1A). 

S phase parameters following Gl block-and-release of 
S. cerevisiae cells 

The progression through S phase of 5. cerevisiae cells 
released from a-factor-mediated Gl arrest was monitored 
by flow cytometry (Figure IB). The relative amount of 
DNA, x, was set to x = 1 in Gl cells. A broad peak of 
fluorescence intensity was observed at all time points 
because of imperfect synchrony and because the PI fluor- 
escence follows a Gaussian dispersion whose standard de- 
viation varies with DNA concentration as a = C*/x, 
where C is a constant (C = 0.24 ± 1.3 x 10" 3 ) (39). As 
the a-factor was washed out, the cells entered S phase 
and the DNA content increased to reach x = 2 (G2/M 
phase). To measure the starting and ending times of S 
phase, we compared the variation in the shape of two 
consecutive flow cytometry profiles by calculating their 
linear correlation coefficient (R) and the associated P 
value (Figure 1C). A coefficient of 1 corresponds to two 
identical curves. The more two curves differ, the lower 
their correlation coefficient. We observed (Figure 1C) 
that R started to decrease between 15 and 30 min 
(R = 0.93, P< 10" 10 ), dropped sharply between 30 and 
45 min (R = 0.4, P = 0.03) and reverted to R = 0.99 
between 60 and 75 min (P < 10~ 10 ). Therefore, cells pro- 
gressed through S phase between 30 and 75 min after 
a-factor release. The number of Gl cells decreases 
sharply between 30 and 45 min (Figure IB), suggesting 
that the rate of S phase entry was maximum over this 
period. To measure the length of S phase, the correlation 
coefficient changes was interpolated by a Gaussian 



are fitted to equation (1) using a simplex algorithm (solid line, 
X 2 =1.03 and /><10~ 4 ). (C) Linear correlation coefficient between 
two consecutive flow cytometry profiles. The linear correlation coeffi- 
cient, R, was calculated between profiles measured at time n— 1 and 
time n. The time n was assigned to the obtained value of R. The 
flow cytometry profile at time = 0 was compared with itself and 
therefore R = 1. The starting time (t aliM = 30 min) and the length of 
S phase (Ts — 42 ± 15 min) were measured by interpolating the 
obtained profile to a Gaussian function. 
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Figure 2. Analysis of S phase parameters in the same cell population as in Figure 1. (A) Fraction of replicated DNA,/ DNA (i), as a function of time. 
Data points were extracted by measuring the amount of incorporated BrdU as described in 'Materials and Methods' section (open circle). The origin 
of time axis was set to the starting time of the S phase, ( sUlr t, which is 30min after Gl release (time = t— « sta rt)- To show that /dna(') follows a 
sigmoidal shape (8). it was fitted to a Boltzmann sigmoidal function (solid black line, x 2 = 1-003 and P < 10~ 4 ). (B) Total rate of DNA synthesis as 
a function of average DNA content. (C) Time-dependent distribution of DNA content (time = t—t alilrl ). Inset: variation of the distributions width as 
a function of time (time = /— 4 tarl ). (D) Probability, n(t), for a cell to be in S phase as a function of time (time = t— r sUlrt ) calculated as described in 
'Material and Methods' section. The black dashed line is a guide for the eye. 




function (we used a Gaussian function because of the U 
shape of the curve and for sake of simplicity). As flow 
cytometry profiles were taken every 15min, the 
maximum uncertainty of the S phase duration in our ex- 
periment is about 15 min. We estimated the mean length of 
S phase to 42 ± 15 min. 

Using the measured / DNA (7) (Figure 2A), we calculated 
the rate of total DNA synthesis, as a function of DNA 
content (Figure 2B) by differentiating ,/dna(0 as a 
function of time. Then, using the method as described in 
'Materials and Method 1 section, we extracted the time de- 
pendent distribution of cellular DNA content, n(x,t), and 
the fraction of cells in S phase, n(t). The dispersion in 
DNA content was initially narrow but increased as 
S-phase progressed then decreased after 35 min (Figure 
2C, inset). The fraction of cells in S phase at time t, n(t), 
corresponded to the ratio between the integral of n(x,i) 
over x and the total number of cells. The resulting n(i) 
profile is shown on Figure 2D. Note that the start time 
of S phase (/ = 0 in Figure 2) was set at 30 min after Gl 
release. 

The consistency of our analysis was further checked 
by showing that we can reproduce the flow cytometry 
profile of an exponentially growing cell population 
using the /dnaM measured by BrdU incorporation 
(Supplementary Section 1). 




Figure 3. Examples of combed DNA fibres from the same cell popu- 
lation as in Figures 1 and 2. BrdU incorporation is revealed in green 
and total DNA in red. Replication extent is 7% (Early, time after Gl 
release = 60 min), 51% (Mid, time = 90 min) and 98% (Late, 
time = 105 min). 



DNA combing and microarray-based I(t) profiles are 
similar 

We used DNA combing (Figure 3) to measure the 
time-dependent, genome-averaged rate of replication 
origin firing, I(t), and density of replication forks, N t {t). 
The normalized rate of firing, I(t)/I max , and the 
normalized fork density, Nj(i)/N max , were calculated as 
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Figure 4. Saccharomyces cerevisiae temporal profile of normalized rate of origin firing and normalized fork density. I(t) and Nj(f) were normalized 
by their respective maximal values. (A) Saccharomyces cerevisiae temporal profile of normalized rate of origin firing from different experimental data 
sets [inverted triangle, DNA combing data from this work; open circle, microarray data from Yabuki et al. (8); open triangle, microarray data from 
McCune et al. (10); open square, microarray data from Raghuraman et al. (7)]. (B) Saccharomyces cerevisiae temporal profile of normalized 
replication fork density. The symbols are the same as in panel (A). (C) open circle, Population-averaged normalized rate of origin activation 
extracted from DNA combing; filled circle, single-cell normalized rate of origin firing obtained by Wiener Fourier deconvolution as described in 
the text. The solid red line is the 5 points smoothed curve. (D) open circle, Population-averaged normalized fork density extracted from DNA 
combing; filled circle, single-cell normalized fork density obtained by Wiener Fourier deconvolution as described in the text. The solid red line is the 
five points smoothed curve. In all panels, the origin of time axis was set to the starting time of the S phase (time = (— ? s tart)- 



described in the Materials and Methods section and 
compared with those extracted from published 
genome-wide replication timing profiles (7,8,10), as 
described previously (16). Although the yeast strains and 
experimental growth conditions differed between experi- 
ments, all the extracted I(t)/I max and N/(t)/N max profiles 
collapsed together without temporal rescaling (Figure 4A 
and B), except for the data collected by Yabuki et al. (8), 
where the growth temperature differs from the one used in 
other studies. As observed previously (16), all the profiles 
increased during the first two-thirds of S phase and 
decreased before the end of S phase. 

In the DNA combing experiments, the rate of origin 
firing was computed fibre by fibre, then mapped on time 
according to the /dna(0 profile shown on Figure 2A (see 
'Materials and Methods' section for details). Therefore, 
the time origin in the plot of these data (Figure 4A and 
B, blue triangles) is the local time at which each fibre 
started replicating. In contrast, in the case of microarray 
(7,8,10) experiments, the rate of origin firing was 
computed globally over the entire genome and mapped 
to the time in S phase reported in these experiments. 
The fact that the combing and microarray measurements 
of /(/) are in such a good agreement strongly suggests that 
the replication process is spatially homogeneous above the 
mean size of combed fibres (164 ± 5kb). In other words, 
the rate of firing over a domain larger than this 



characteristic length predominantly depends on its extent 
of replication rather than on its absolute starting time. 

The consistency of our data was further checked by 
showing that using a formula that links the rate of 
origin firing to the fraction of replicated DNA (40), we 
can reproduce the /dna(0 curve measured by BrdU in- 
corporation from the I(t) we extracted from microarray 
or combing data (Supplementary Section 2). 

Deconvoluting a single-cell I(t) profile from 
population data 

An important issue in interpreting the population- 
averaged curves shown in Figure 4A and B is that all 
cells do not enter and exit S phase at the same time 
(Figure 2D). This could blur the temporal profile of 
origin firing. However, as the replication process is 
spatially homogeneous, the observed rate of origin 
firing at time t, is the result of origin firing in all 
cells that are in the S phase for a time r < t. In other 
words 1(f) — fl I s (x)n(t — r)dr, where l s (t) represents the 
rate of origin firing in a single cell (assumed to be identical 
among cells) and n (t—r) is the fraction of cells in S phase 
for a time r (r < t). As the kinetic of BrdU uptake by the 
cell is consistent with the flow cytometry profile of an ex- 
ponentially growing cell population (Supplementary 
data), the extracted n(t) can be used to remove the 
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Normalised fork density Normalised fork density 

Figure 5. Normalized rate of origin firing versus normalized fork density extracted from DNA combing data in cell populations (A) and in single 
cells (B). The open black squares (open square) and the open grey circles (open circle) represent the increasing and decreasing part of 1(f), respect- 
ively. The solid black and dashed grey lines are the spline smoothed curves in panel A, and the five points smoothed curves in panel (B). 



contribution of temporal dispersion from the combing 
data using Wiener Fourier deconvolution method. 

The raw data and the resulting single-cell I s (t)/I max and 
N s f(t)/N ma>i profiles are shown on Figure 4C and D, re- 
spectively. The raw data and deconvoluted profiles are in 
very good agreement. The deconvoluted profiles decrease 
at the end of S phase, although less markedly than in the 
raw data. Therefore, although the apparent drop in origin 
firing at the end of S phase is somewhat enhanced by the 
asynchrony of the cell population, it is still observed at 
the single-cell level. 

The shape similarity of the I(t) and Nj(t) profiles, which 
both peak at 30min in S phase, demonstrate a correlation 
between fork density and rate of origin firing. This is 
further exemplified by plotting the normalized rate of 
origin firing against the normalized replication fork 
density for the DNA combing experiment (Figure 5A). 
Removing temporal dispersion as described above strik- 
ingly improved the observed correlation (Figure 5B). 

Modelling l(t) profiles in wild-type and clb5A 
mutant strains 

To push this analysis further, we extracted and compared 
the rate of origin firing and fork density profiles for the 
wild-type strain and the clb5A mutant studied by McCune 
et al. (10) (Figure 6A and B). The rate of origin firing and 
the fork density in the clb5A mutant increased and 
decreased more slowly than in the wild-type and their 
maxima were 40% lower. Interestingly, the deletion of 
CLB5 drastically reduced the proportionality between 
rate of origin firing and fork density during late S phase 
but not during early S phase (Figure 6C). The crossover 
between the two regimes occurred when /(/) reached its 
maximal value. We previously demonstrated that a 
minimal model of stochastic origin firing could quantita- 
tively describe the time-dependent changes in rate of 
origin firing and fork density observed in a wide range 
of eukaryotes including the yeast S. cerevisiae (16). In 
this model, origin firing is governed by the encounter of 



a recyclable replication-fork factor with a potential origin 
to produce an firing event with probability P(t). The 
amount of this factor increases during S phase at rate 
J and the probability P(t) increases with fork density ac- 
cording to a self-limiting law P = Pq+P\[\ — e~ NbA ^ c ], 
where A^f or k is the time-dependent number of forks and 
N c the maximum number of forks. A simulated 
normalized rate of origin firing (average of 1000 simula- 
tions) with J = 5 factors per calculation round and 
P=3x \Q- 4 +2 x 10- 4 x [1 - e-AWWoO] is shown on 
Figure 6D, grey dashed line. The obtained I{t) profile is 
similar to the experimental profile for the wild-type strain 
(Figure 6 A, grey dashed line). We then explored how this 
model needs to be modified to account for the experimen- 
tal I(t) profile of the clb5A mutant. This mutant shows a 
longer S phase associated with the following changes in 
I(t) profile: (i) the ascending part stops before reaching the 
maximum value of the wild-type; (ii) the peak is replaced 
by a plateau; (hi) the correlation between origin firing and 
fork density is abolished once the plateau is reached 
(Figure 6A and C). Therefore, we modified the model by 
allowing the rate / and the dependence of the firing prob- 
ability P(t) on the fork density to decay exponentially: 
/= 5 x e -W28 and p = 3 x io- 4 +2 x 10" 4 x e ~'™ /28 x 
[1 — e~ Nfotk / U0 °], where t simu is the number of calcula- 
tion rounds (Figure 6D). The obtained profile strikingly 
reproduced the observed behaviour of the rate of ori- 
gin firing in the clb5A mutant. Importantly, the I(t) 
profile of the clb5A mutant could solely be reproduced 
if both / and P 1 were forced to decay. Therefore, 
although our model is much simpler than the complex- 
ity of the underlying biochemical circuitry that controls 
origin firing, it can account for the behaviour of both 
the wild-type and clb5A strains. It is possible that 
errors on detection of replication origins and their 
firing time might differently affect the wild-type and 
clb5A I(t) profiles. Nevertheless, we believe that such 
errors are not ample enough to significantly affect our 
conclusions. 
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Normalised fork density Time (round of calculation) 

Figure 6. Modelling wild-type (dashed grey) and clb5A (solid black) rate of origin firing and fork density. (A and B) Normalized rate of origin firing 
(A) and fork density (B) from wild-type (open circle) and clb5A (open square) strains, extracted from McCune et al. (10). All data sets were 
normalized to the maximal value obtained for the wild-type strain. (C) Plot of the normalized rate of origin firing versus normalized fork density in 
the wild-type and clb5A strains. The increasing and decreasing parts of I(t) are indicated by arrows. (D) Simulated normalized rate of origin 
activation as a function of time. The dashed grey line represents the unmodified model (wild-type [WT]) and the solid black line the modified 
model (c!b5A). 



DISCUSSION 

It has been debated whether eukaryotic DNA replication 
origins fire independently of each other or whether the 
firing of late origins is tied to the progression of earlier 
assembled forks. Two important kinetic parameters of S 
phase that are relevant to this question are the 
genome-averaged, time-dependent rate of origin firing, 
I(t) and replication fork density, N/f). Recent studies 
suggest that I{t) has a universal shape among eukaryotes. 
A model that explicitly ties I(t) to Nj(f) has been proposed 
to explain this shape. However, some measurements of 
and N/(t) were based on population averages, which 
may introduce errors due to imperfect cell synchrony 
and cell-to-cell variability. Here, we have computed 
population-averaged I(t) and N/(t) profiles from micro- 
array data and from DNA combing data, and we have 
extracted single-cell I s (t) and N s j(t) profiles, by 
deconvoluting from the combing data the fraction of 
cells n(t) in S phase at time / in the same cell population. 
We find that all three I{t) profiles and all three N f (i) 
profiles are in good agreement. 

The similarity of the profiles computed from the 
combing or from the microarray data indicates that 
errors specifically associated with interpretation of each 
dataset are unlikely to strongly affect the shape of these 
profiles. 

The combing /(/) and Nf{t) profiles were obtained by 
computing the firing rate /(/dna) and fork density 
N/frmjO as a function of replication fraction / DNA and 
by converting to time using the /dna(0 curve determined 
by analysis of BrdU incorporation in the same cell 



population. This strategy implicitly assumes that I(t) and 
Nj(t) are spatially homogeneous. As combing is done on 
short DNA fragments (165kb on average) and is consist- 
ent with genome-averaged results obtained from micro- 
array data, we conclude that it is unlikely that /(/) and 
Nj(i) widely vary over the genome, and that the spatial 
homogeneity assumption is valid. This does not argue 
against the existence of domains of the genome that are 
preferentially replicated early or late during the S phase 
(10). However, our results suggest that in a general 
manner the rate of origin firing and the density of forks 
in a particular locus of the genome strongly depends on 
the local degree of replication of this locus. 

By deconvoluting «(/) from I(t) and Nf(t) the population 
distribution contribution is removed from the profiles. 
The similarity of the raw and deconvoluted profiles 
argues that they are an intrinsic property of the DNA 
replication process and are only mildly affected by the 
imperfect synchrony and cell-to-cell variability of the cell 
population. The deconvoluted profiles do decrease at the 
end of S phase, although less markedly than in the raw 
data. This indicates that population averaging exagger- 
ates, but does not artefactually create, the drop in I(t) 
and N/(t) at the end of S phase. 

We demonstrated a linear correlation between and 
N/{t). Although it is trivial to have I(t) ~N,(t) at the be- 
ginning of S phase, as each initiation creates two forks, it 
is not trivial to have such a correlation at the end of S 
phase, because termination reduces N f {t) and a priori does 
not affect initiation. In principle, it is possible that both 
Nf(i) and I{t) are controlled by a third function. But as 
N/j) can only be affected by initiation and termination, 
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this third function would have to regulate termination 
by acting on fork velocity, which is difficult to reconcile 
with the reported uniformity of fork progression through 
the genome (35). Therefore, we propose that N/f) 
controls I(t). 

This proposal is based on late S-phase DNA combing 
data. A potential concern with DNA combing is that an 
origin that fired close to a fork will generate a small eye 
that can quickly merge with the neighbouring eye and 
escape detection. This bias will increase through S phase 
as there are more and more eyes to merge with. To address 
this concern, we ran a set of simulations that show that the 
linear relation between I{t) and Nj(t) is independent of the 
resolution with which we can detect replication eyes and is 
only observed when Nf(t) controls I(t) (Supplementary 
Section 3). 

Interestingly, the linear correlation between Nj{t) and 
I(t) is drastically reduced in late S phase, but not in 
early S phase, in the clb5A mutant. The correlation is 
lost once I(t) reaches its maximal value. While in the 
wild-type /(/) reaches a maximum, in the mutant I{t) 
reaches a plateau, at 60% of the wild-type maximum, 
that reflects a limitation of initiation over a longer 
period of time. It has already been reported (27) that in 
the clb5A mutant late origins fire less efficiently, thus, 
inducing the lengthening of S phase. 

We previously proposed a minimal model of stochastic 
origin firing that could quantitatively describe the Nj(t) 
and I(t) profiles observed in a wide range of eukaryotes. 
In this model, I(t) is governed by the encounter of a re- 
cyclable replication-fork factor with a potential origin to 
produce an firing event with probability P(t). The amount 
of this factor increases during S phase at rate J and the 
probability P(t) increases with Nj(t) in a self-limiting 
manner. We found that the difference in wild-type and 
clb5A mutant N/(t) and /(/) profiles can be reproduced if 
both J and the coupling factor Pj between N/J) and P(t) 
decrease exponentially as the simulation progresses. 

The need to make J decay is straightforward to inter- 
pret. During a normal S phase, CDK activated by Clb5p 
or by Clb6p stimulates origin firing. The amount of Clb5p 
protein is stable during S phase but Clb6p is gradually 
degraded at the beginning of S phase (24,25) so that in 
the absence of Clb5p the S phase is lengthened due to a 
decay in CDK activity. CDK activity is required to 
promote the assembly of replication forks from multiple 
protein factors including some that are limiting for origin 
firing. 

The need to make P\ decay exponentially is justified by 
the looser dependency of /(/) on N f (t) in the clb5A mutant 
(Figure 6C). Interestingly, late replicating regions that are 
affected in the clbSA mutant are affected as well in 
mutants deficient in the chromatin remodelling factors 
Isw2 and Ino80, which are enriched in replicating 
regions in a similar manner to Pol a (28). Therefore, we 
speculate that fork-associated chromatin remodellers may 
be specifically required to promote origin firing. In other 
words, fork-stimulated origin firing would become pre- 
dominant in late-replicating regions and would be ineffi- 
cient in the clb5A mutant due to the decreased affinity of 
limiting initiation factors for potential origins in the 



presence of lower CDK activity. Therefore, our results 
suggest that any modifications in chromatin state (struc- 
ture or chemical state) would influence the dynamics of 
DNA replication and particularly the position and firing 
time of replication origins. This hypothesis is consistent 
with a direct link between epigenetic modification of chro- 
matin and replication timing. 
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Supplementary Data, Supplementary Figures 1-4, 
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